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The experimental realization of lattices with Chem bands in ultracold-atom and photonic systems has moti¬ 
vated the study of time-dependent phenomena, such as spatial propagation, in lattices with nontrivial topology. 
We study the dynamics of gaussian wavepackets on the Haldane honeycomb Chern-band lattice model, in the 
presence of a harmonic trap. We focus on the transverse response to a force, which is due partly to the Berry 
curvature and partly to the transverse component of the energy band curvature. We evaluate the accuracy of a 
semiclassical description, which treats the wavepacket as a point particle in both real and momentum space, in 
reproducing the motion of a realistic wavepacket with finite extent. We find that, in order to accurately capture 
the wavepacket dynamics, the extent of the wavepacket in momentum space needs to be taken into account: The 
dynamics is sensitive to the interplay of band dispersion and Berry curvature over the finite region of momen¬ 
tum (reciprocal) space where the wavepacket has support. Moreover, if the wavepacket is prepared with a finite 
initial momentum, the semiclassical analysis reproduces its motion as long as it has a large overlap with the 
eigenstates of a single band. The semiclassical description generally improves with increasing real-space size 
of the wavepacket, as long as the external conditions (e.g., external force) remain uniform throughout the spatial 
extent of the wavepacket. 


I. INTRODUCTION 

The presence of non-trivial topology in the energy bands 
of lattice models [1, 2] has generated intense interest over 
the past decade. When an energy band possessing a nonzero 
Chern number is filled with fermions, one obtains a Chern 
insulator, realizing a quantized Hall effect without external 
magnetic fields. The Chern number is a topological invariant, 
defined as the flux of the Berry curvature over the Brillouin 
zone. The Berry curvature acts like a momentum space ana¬ 
logue of the magnetic field [3]. One of the first examples of 
a Chern band model was introduced in the seminal work by 
Haldane [4], where time reversal symmetry is explicitly bro¬ 
ken in a honeycomb lattice model of free fermions by complex 
hoppings between next-nearest-neighbor lattice sites. 

Chern band physics has been discussed in numerous dif¬ 
ferent scenarios. Originally, the focus was on electronic sys¬ 
tems [5-9]. More recently, lattices with nontrivial Berry 
curvature have been experimentally realized using ultracold 
atoms trapped in optical lattices [10-13]. The spatial geome¬ 
tries in these setups are more flexible than in traditional solid- 
state situations. These experimental developments have thus 
motivated the study of effects of lattice and trap geometries on 
topological states [14-19]. Cold atoms also provide an excel¬ 
lent platform for observing and analyzing non-equilibrium dy¬ 
namics of Chem bands. Therefore, there is increasing interest 
in developing protocols to characterize the topological nature 
of Chern bands from non-equilibrium behavior [20-30] . For 
instance, the quench dynamics of Chern insulators can signal 
their non-trivial edge states [26] and the semiclassical trajec¬ 
tory of wavepackets can be related to the Chern number [22]. 

Topological photonics has proved itself to be another im¬ 
portant context for Chern band physics [31-40]. Theoretical 
studies and proposals [32-36] have also been supplemented 
by experimental observation of the topological edge states in 
photonic systems [37-40]. In addition, topological magnons 
have been proposed as a platform to study transport influenced 
by Berry curvature in the context of the thermal Hall effect 
[41-43]. 


For ultracold atoms (especially bosons) and for photonic 
systems, a dynamical situation where the atoms/photons form 
a spatially localized and evolving wavepacket is more natu¬ 
ral than a static situation in which a band is exactly filled. 
Thus, the recent experimental developments call for a system¬ 
atic understanding of the non-equilibrium dynamics of real¬ 
istic wavepackets on Chern band models. A particularly im¬ 
portant theme is the response of a localized wavepacket to an 
applied force (potential gradient) [13, 22, 27, 29, 44]. One ex¬ 
pects Bloch oscillations in the direction of the applied force. 
In addition, there is also motion perpendicular to the direction 
of the force, with contributions due to the topological Berry 
curvature and due to the band dispersion. 

A natural framework to describe the motion of wavepackets 
is to use semiclassical equations of motion. For systems with 
Berry curvature, the semiclassical equations were introduced 
in the context of crystals with magnetic Bloch bands [45-50], 
anomalous Hall responses [50-52], and optical lattices with 
spin-orbit coupling [53]. They have been used to study the ef¬ 
fect of the Berry curvature on wavepacket trajectories [29, 54] 
and on collective modes [55, 56] in ultracold gases. For in¬ 
stance, Ref. [29] outlines a procedure to isolate the contribu¬ 
tion of the topological Berry curvature by separately evolving 
the system under opposite potential gradients and then appro¬ 
priately summing up the responses. 

Semiclassical approaches typically rely on the approxima¬ 
tion of assigning a sharply defined position and momentum 
simultaneously to a quantum state. Based on this assump¬ 
tion one can solve the set of coupled differential equations 
for position and momenta and hence obtain sharply defined 
trajectories of the particle in both real and momentum space. 
However, this assumption is a priori not valid in realistic situ¬ 
ations where the wavepackets have a finite spread in both real 
and momentum space, a scenario expected in typical ultracold 
bosonic and photonic experiments. 

In this work, we study the effect of size and initial momen¬ 
tum of a wavepacket placed off-center in a harmonic trap in 
Haldane’s honeycomb model, and evaluate the ability of the 
semiclassical approach to describe these effects. We focus on 
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time evolutions up to moderate timescales such that the dis¬ 
placement of the wavepacket from its initial position is not 
large compared to the wavepacket size and much smaller than 
the distance to the center of the trap. In this regime, the trap 
plays a role similar to a constant force. We quantify the trans¬ 
verse motion of the wavepacket using the angular velocity 0 
with respect to the center of the trap. 

We find that the standard point-particle semiclassical ap¬ 
proach captures some qualitative features of the dynamics, but 
is generally insufficient to quantitatively describe the actual 
real time dynamics. Therefore, we reformulate the semiclas¬ 
sical description to take into account the finite spread of the 
wavepacket in momentum space. The exact evolution of 0 
is compared in detail to predictions from the extended semi¬ 
classics and from the standard point-particle semiclassics. In 
momentum space, the wavepacket moves at a constant rate in 
the direction of the force, leading to Bloch oscillations. As it 
crosses different regions of the Brillouin zone, the local Berry 
curvature and band curvature determine 9{t). For spatially 
localized wavepackets, the extent in momentum space is fi¬ 
nite. The “extended semiclassics” procedure incorporates the 
variations of band dispersion and Berry curvature in this ex¬ 
tended region of momentum space. We find that, as long as 
the physics is dominated by one band, this procedure repro¬ 
duces the full dynamics well. This shows that the basic idea 
of semiclassics (simultaneously assigning both position and 
momentum to a quantum particle) can incorporate aspects of 
the full quantum dynamics to an extent beyond what is known 
from the usual point-particle treatment. 

One might intuitively expect that semiclassical descrip¬ 
tions should work better for spatially large wavepackets, since 
these correspond to smaller regions in momentum space. We 
show that this is generally true, but that semiclassics still de¬ 
scribes the dynamics of rather small wavepackets, especially 
if momentum-space extent is included. In addition, by con¬ 
sidering a tight trap, we show an example of possible experi¬ 
mental relevance where larger real-space sizes can render the 
semiclassical description less inaccurate, due to an inhomo¬ 
geneity of the force within the spatial support region of the 
wavepacket. 

We also demonstrate the effect of initializing the 
wavepacket with a finite momentum. In addition to zero mo¬ 
mentum (T point) we start with the packet at one of the K 
points and one of the M points of the Brillouin zone. We 
show that imparting momenta to a gaussian wavepacket us¬ 
ing an e®'^ ''-like factor can cause the wavefunction to have 
significant occupancy in the upper band, including extreme 
cases where it is almost completely transferred to the upper 
band. As long as one of the bands dominates, the semiclas¬ 
sical description works well when using the properties (band 
dispersion and Berry curvature) of the band where the state 
has most of its weight. The single-band semiclassical pro¬ 
cedure is naturally insufficient when multiple bands are sig¬ 
nificantly occupied: Features like interference oscillations are 
not captured by an incoherent averaging of contributions from 
different bands. 

The rest of the paper is organized as follows: in Sec.II we 
describe the model, the geometries and the simulations of time 



FIG. 1. (Color online) (a) Honeycomb lattice defining the Haldane 
model (1) with black and grey circles denoting the two sublattices. 
The solid lines represent the real nearest-neighbor hoppings whereas 
the dotted lines the complex next-nearest-neighbor hoppings. The 
arrows indicate the direction in which the particle hops to pick up 
a phase <)>. (b) The two band energy spectrum with parameters 

Ji = — 1, J 2 = —1/4, and (j) = 0.497r. (c) The Berry curvature 
in the lower band for the same parameter values shown as a color 
map in the Brillouin zone. The high symmetry points are marked. 
The boundaries of the Brillouin zone are formed by the reciprocal 
lattice vectors Gi = (0, 47r/3) and G 2 = (27r/\/3, —27r/3). 


evolution. In Sec. Ill we discuss the semiclassical framework 
for the dynamics of the wavepacket, followed by a compari¬ 
son of the results from simulations and semiclassics in Sec.IV. 
Finally, we comment on the dynamics in a tight harmonic trap 
in Sec.V and provide some discussion and context in Sec.VI. 


II. SIMULATIONS OF WAVEPACKET DYNAMICS ON 
THE HALDANE MODEL 

A. Model Hamiltonian 

Haldane’s model [4] is a tight-binding Hamiltonian on a 
honeycomb lattice: 

+ J 2 SJSj-f h.c (1) 

(ij> 

where the (I, j) denotes the nearest neighbors, {{i,j)) denotes 
the next nearest neighbors and S] (bj) is the creation (annihi¬ 
lation) operator at site i. We will formally consider the dy¬ 
namics of a single particle, so fij, bj may be thought of as 
either fermionic or bosonic operators. When particles hop be¬ 
tween next nearest neighbors they pick up a phase = (p 
if they hop in the direction of the arrow shown in Fig. 1(a) 
and —p if they hop in the opposite direction. The energy 
spectrum is gapped if p ^ mr and is particle hole symmet¬ 
ric for p = tt/2. Throughout this paper we use the parameters 
Ji = —1, J 2 = —1/4, and p = 0.497r, for which the energy 
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spectrum is depicted in Fig. 1(b). The parameters are cho¬ 
sen so that the two energy bands are quite similar and none of 
them are excessively flat (since band flatness can introduce ad¬ 
ditional peculiarities in the dynamics), and so as to avoid ex¬ 
act particle-hole symmetry, since we are interested in generic 
rather than fine-tuned effects. We expect the physics described 
in this paper (explicitly for these parameters) to be exemplary 
for a wide region in parameter space. 

For these parameters the bands have a finite Chern number; 
— 1 for the lower band and -|-1 for the upper band. The distri¬ 
bution of the Berry curvature of the lower band in the Brillouin 
zone is shown in Fig. 1(c). The upper band has approximately 
opposite Berry curvature, i.e., positive instead of negative val¬ 
ues. We provide some details of the topological properties and 
conventions in Appendix A. Fig. 1(c) also shows the high- 
symmetry points. In addition to zero-momentum (F point), 
these are three inequivalent M points and two inequivalent K 
points. We choose the boundaries of the Brillouin zone to be a 
parallelogram. In the literature this is equivalently often cho¬ 
sen to be of hexagon shape; the solid lines inside the Brillouin 
zone show the boundaries for such choice. 

We set h = I, measuring time in units of h/Ji and energy 
in units of Ji. Space and momentum are in units of lattice 
spacing (set to unity) and inverse lattice spacing, respectively 
and geometric angles are measured in radians. 


B. Construction of wavepackets 

We are interested in the dynamics of a wavepacket of finite 
real-space extent in Haldane’s honeycomb model in the pres¬ 
ence of a harmonic trap. The initial wavepacket is prepared 
with gaussian shape. For zero initial momentum, 

\^(t = 0)) = |;). (2) 

Here \l) denotes a single-particle state, with the particle com¬ 
pletely localized at a site indexed by I and JV is a normal¬ 
ization factor, and r; = {xi,yi) denotes the Euclidean posi¬ 
tion of site 1. The coefficient c;(f) denotes the amplitude of 
the wavefunction at time t at site I, and a is the width of the 
wavepacket. We will also use wavepackets with nonzero ini¬ 
tial momenta. A finite momentum is obtained by multiplying 
the coefficients with a site-dependent phase factor; 


well-resolved in order to obtain the energy-momentum disper¬ 
sion. In our case, we supply the wavepacket with a momentum 
using the factor but do not specify energy. This can be 

thought of as a Bragg pulse with sharp momentum resolution 
but poor or non-existent energy resolution. This allows us to 
explore various occupancies of the two energy bands. 

The momentum zero wavepacket (2) turns out to predomi¬ 
nantly have overlap with eigenstates of Hum at the bottom of 
the spectrum, in the lower band, as long as a is not too small. 
This is generally true in simple lattice models with negative 
hopping constants. In a complicated model like Hum, this is 
not a priori obvious, but is the case for the parameters we are 
using. 

Boosting the wavepacket in momentum space as in Eq. (3) 
can result in the wavepacket having support on both the bands 
of "Hhm- This is exemplified in Eig. 2 through the overlap of 
the wavepacket with the eigenstates of 'Hhm- We denote the 
overlap of the initial state, = 0)) with an eigenstate of 
'Hhm, l^a) with eigenvalue as Oa = — 0)|uq,)P. 

Eig 2(a) shows a plot of Oa against Ea for four out of the six 
high-symmetry momentum points being the initial momen¬ 
tum of the wavepacket. It can be seen that for ko = K the 
wavepacket has support on both bands. The wavepacket corre¬ 
sponding to ko = Ml is also shifted higher in energy, though 
it has overlaps primarily with the states of the lower band. On 
the contrary, for ko = M 2 the weight shifts almost entirely to 
the upper band. Such a drastic difference of behavior between 
the Ml and M 2 points may seem unexpected because they 
are related by symmetry. However, the eigenfunction struc¬ 
tures are of course inequivalent, so the overlap distributions 
after a momentum boost cannot be expected to be similar. 

We quantify the weight of the wavepacket on the lower 
band, W- by taking the sum of the overlaps of the wavepacket 
with the eigenstates of the lower band. Mathematically, 

N/2 

IE_ = ^0„; W+ = l-W., (4) 

a—1 

where is the weight on the upper band, and N is the 
number of sites in the lattice and hence the number of single¬ 
particle eigenstates. The color map in Eig. 2(b) shows the 
magnitude of W- for a wavepacket with a given initial mo¬ 
mentum in the Brillouin zone. 


C. The trap 


\tp{t = 0; ko)) = \l) 




\ri - Tc 

2a2 


iko • r; 


10 - ( 3 ) 


The evolution in time is carried out with the Hamiltonian 
T~L = 'HnM + ’Htrap = 'Hhm + ^ (5) 

i 


In cold atom experiments, a wavepacket or atomic cloud can 
be boosted in momentum space in this way using a ‘Bragg 
pulse’; this is commonly used to determine the excitation 
spectrum of cold atom systems using Bragg spectroscopy (see, 
e.g.. Ref. [57] for a description of the experimental technique). 
In experimental spectroscopy, both energy and momentum are 


where 

E(0 = ^|ri-roP, (6) 

with To = (xo,yo) denoting the center of the harmonic trap 
and Vo controlling its strength. 
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FIG. 2. (Color online) (a) The overlaps Oa of the initial wavepacket 
at different initial momenta ko with the eigenstates of HmA (with 
periodic boundary conditions) are plotted against the energy eigen¬ 
values Ea. (b) The weight of the wavepacket on the lower band 
{W-) for each initial momentum in the Brillouin zone. 


The force exerted by the trap is along the inward radial di¬ 
rection, so that one expects Bloch oscillations in this direction. 
We are particularly interested in the transverse response, and 
hence the angular velocity of the wavepacket around the cen¬ 
ter of the trap is a natural observable to study. The angular 
variable 0 at position r = (x, y) is given hy 6 = tan“^ x-xl ■ 
Its average as a function of time is calculated using 


{0){t) = tan ^ 


{y){t)-yo 

{x){t) - Xo’ 


( 7 ) 


from the average x and y for the time dependent wavefunction 


{x)it) = '^\ci{t)\‘^xi, {y)it) ='^\ci{t)\‘^yi. (8) 

i i 


D. Simulations of wavepacket dynamics 

In this work we present results for the exact dynamics of 
wavepackets placed off-center in the trap, and compare with 
semiclassic predictions. 

We present simulations mostly for a relatively weak trap 
(Vo = 0.002), where the wavepacket width a is much smaller 
than the distance (« 164.5) to the trap center rg. The trap 
potential gradient (i.e., the force) does not vary too much over 
the extent of the wavepacket. In Sec. V we also present results 
for a tighter trap (Vb = 0.02), where the initial distance of 
the wavepacket to the trap center is 10 times smaller so as to 
have the same force at the center of the wavepacket. The trap 
curvature is more significant in this case. 

In Fig. 3(b) we show the real-space evolution of the 
wavepacket in the weak trap, for two different initial mo¬ 
menta. We focus on dynamics up to f « 20. The motion 
of the wavepacket on this timescale is not large compared to 
its width. The force acts in the radial (positive x) direction. 
A transverse response, perpendicular to the force, is clearly 
visible; we analyze this quantitatively through the time depen¬ 
dence of {9). In the next sections we provide a thorough com¬ 
parison of numerically exact results for {9) obtained through 


direct simulation (which we refer to as (0) exact) with predic¬ 
tions from the semiclassical formalism, to be defined below. 

In addition to the transverse response, there are also Bloch 
oscillations in the radial direction; this is not obvious in the 
real-space snapshots but is more evident in momentum space. 
The motion of the wavepacket in momentum space is obtained 
by taking a lattice Fourier transform of the coefficients C; (t) 
at each instant of time to obtain the occupancies of each mo¬ 
mentum, denoted by CK(f). This motion can be visualized by 
plotting the coefficients |ck(f)P over the Brillouin zone at dif¬ 
ferent instants of time, as done in Fig. 3(c). The wavepacket 
moves through the Brillouin zone at constant rate in the di¬ 
rection of the force. Due to the periodicity of the Brillouin 
zone, each time the wavepacket exits through the right or top 
boundary it re-enters through the left or bottom boundary. As 
a visual aid, we show schematically in Fig. 4 the trajectories of 
the wavepacket centers starting from the three high-symmetry 
points. 


III. SEMICLASSICAL DYNAMICS 


In this section, we set up the semiclassical framework to 
calculate the time evolution of 0(t). We first formulate the 
basic ‘point-particle’ approach, under the standard assump¬ 
tion of simultaneously well-defined position and momentum. 
We then formulate an extension where the structure of the 
wavepacket in momentum space is taken into account. 

In the most basic semiclassical approach, the structure of 
the wavepacket in both real and momentum space are ignored, 
so that the wavepacket or particle is described by a sharply de¬ 
fined position and momentum (r, k). In addition, it is also as¬ 
sumed that the wavepacket dynamics is completely governed 
by a single energy band. We will first write down the semi¬ 
classical equations with the assumption that the wavepacket 
has support only on the lower band. The semiclassical equa¬ 
tions of motion are [50] 

nV nli 



(9b) 


Here E- (k) is the energy dispersion, and r2_ (k) is the Berry 
curvature of the lower band. From the second term in Eq. (9a), 
we see that the Berry curvature induces a velocity perpendic¬ 
ular to the direction of the external force, which leads to the 
transverse motion of the wavepacket. 

We now specialize to the geometry we are using, with a 
trapping potential centered at {xo,yo)- Using (6), the external 
force is given by 


F(r) = -VrU(r) = -Vo{{x - xo)x + {y - yo)y). (10) 


So, the semiclassical equations form a set of four coupled dif¬ 
ferential equations 

t;_,,(k) = ^ = ^L_^ + Voiy-yo)niCk), (11a) 
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FIG. 3. (Color online) (a) Geometry of the lattice and the trap. The surface represents the potential of the harmonic trap and the colors show 
the occupancies of the wavepacket at f = 0 on the real space lattice. The lower figure shows a zoom near the wavepacket. (b) Snap shots 
of the time evolution of the real space occupancies (|c;(f)p). The red circle shows the center of the mass of the wavepacket. The solid lines 
are equipotential contours of the trap potential, (c) Occupancies of each momentum mode in the Brillouin zone at different times (|ck(f)P). 
The trap strength for all figures is Vq = 0.002. The initial wavepacket has a gaussian width = 50, initial momentum ko = M 2 , and 
It’D — fcl = 164.5 in units of the lattice constant. 



FIG. 4. (Color online) Trajectories of the wavepacket in momentum 
space for three initial momenta showing Bloch oscillations. The or¬ 
ange circles represent the trajectory of the center of the wavepacket 
in momentum space. The blue larger circles denote the position of 
the wavepacket center in momentum space at four instants of time 
{t = 0, 5, 10, and 15). 


This set of equations can be solved explicitly to trace out the 
trajectory in time of a particle in real as well as momentum 
space. We label as (0)_^pp.sc the angular velocity correspond¬ 
ing to the real space trajectories calculated in this way. (The 
subscript ‘pp-sc’ stands for ‘point-particle semiclassics’ and 
the — sign denotes that the lower band properties have been 
used). A similar calculation can be done with the character¬ 
istics of the upper band (£)_(. (k) and G^(k)) and the angular 
velocity so calculated is denoted by (0)_|_ pp.sc. 

As observed previously, wavepackets can have support on 
both bands. One reasonable procedure would be to use the 
(0(f))_,pp.sc or (0(f)) +, pp-sc curve, depending on whether the 
lower or upper band has more occupancy. We follow a some¬ 
what more refined procedure by taking the weighted average 
of the two according to the weights of the initial packet 
on the two bands. Hence the angular velocity calculated from 
the point-particle semiclassics is defined as 


(0)pp-sc = H^-(0)-.pp-sc + kk+(0) + ,pp-sc (12) 


u_,y(k) = ^ = - V))(a: - a;o)0^(k), (11b) 


=-Vq[x-xq), (11c) 

at 

=-Vo{y - yo). (lid) 

at 


A key assumption above is that the wavepacket can be 
treated like a point particle in both real and momentum space 
simultaneously, hence neglecting the quantum nature of the 
wavepacket. However in realistic quantum experiments and 
simulations, where the wavepacket is of finite extent, the va¬ 
lidity of this assumption is not a priori clear. We now ex¬ 
tend this formalism to take into account the finite spread of 
the wavepacket in momentum space. From the geometric def- 
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inition 9 = tan ^ H—M^ ^ve obtain 

X — Xq ’ 

■ {x - Xo)vy - {y - yQ)v^ 

{x - xof + (y - yoY ’ ^ ^ 

By using the expressions of t^±,x(k) and y(k) obtained 
from Eq. (11a) and Eq. (11b), we define the functions 9± (k) in 
the Brillouin zone. Their typical prohles are shown in Eig. 5, 
with parameter values corresponding to the initial position 
used in Eig. 3. 

With our parameters, we have f?-(k) « —i?_|_(k). Also, 
we always have fl+(k) = —fl_(k). Hence we get 0_(k) « 
—0+(k). In other words the prohles shown in Eig. 5 for the 
two bands, 0±(k), are nearly but not exactly negative of each 
other. 


ijkl 



0 7r/2 TT 





0.008 


0.004 


0.000 


-0.004 


-0.008 


0 7r/2 


FIG. 5. (Color online) Profiles of 0-(k) and ^+(k) in the Brillouin 
zone as calculated from combining the semiclassical equations of 
motion Eq. (9a) and the kinematic relation in Eq. (13). The coordi¬ 
nates of the center of the mass of wavepacket in real space relative to 

the center of the trap are given by I—= —164.5 and jz—r/o = 0.5, 

and the trap strength is Vo — 0.002. Here 0-(k) and S+(k) have 
nearly but not exactly equal and opposite values, 0-(k) « —0+(k). 


These profiles of 6± (k) can be used to calculate the evo¬ 
lution of the angular velocity in time by taking a weighted 
average of 0±(k), the weights being the occupancies of the 
wavefunction in momentum space (|ck(f)P) multiplied by the 
weights in each band iW±) defined in (4). We denote the an¬ 
gular velocity calculated this way as (0)wp-sc(f), the ‘wp’ as a 
reminder that the wavepacket structure is taken into account. 
Thus 


(0(f))wp-sc = VE-(0(f))-.wp-sc + kE+(0)+.wp-sc 

= H^_ ^ |5k(f)p0-(k) + kE+ 5] |£k(f)p04-(k) 

kGBZ kGBZ 


^ |ck(f)p[vE-0_(k) + IE+0+(k) 


kGBZ 


(14) 


This procedure assumes that the force does not change along 
the spatial extent of the wavepacket; the spread of the 


wavepacket in momentum space is taken into account while 
a point-particle description is used in real space. Therefore, 
this description will break down when the force varies signifi¬ 
cantly within the real-space support region of the wavepacket 
(Sec. V). In addition, note that this is not a computationally 
advantageous approximation for the time evolution, since we 
are anyway solving the full problem in order to obtain the 
Eourier transform Ck(f). Our purpose here is to investigate 
whether (and how much) taking the momentum-space spread 
into account improves the semiclassical description. 

In this work we focus on parameter regimes such that the 
wavepacket does not have large displacements in real space 
within the time scales f < 20 of interest (Eig. 3). There¬ 
fore, we make a further simplifying assumption and take 9 
as position independent, setting r to be initial position of the 
wavepacket at f = 0, and use the resulting distribution of 9 
to calculate the average. In the following two sections we test 
how well this procedure describes the angular motion of the 
wavepacket. 


IV. COMPARISON OF SEMICLASSICAL PREDICTIONS 
WITH EXACT DYNAMICS 

In this section, we compare the angular velocities of the 
wavepacket obtained from the exact simulations, ((0)exact) to 
those obtained from the semiclassical calculations ((0)pp-sc 
and (0)wp-sc) and discuss the regimes of validity of the semi¬ 
classical framework. 

In Eig. 6 we plot the angular velocities as a function of time 
for the setup corresponding to that shown in Eig. 3(a) for two 
different sizes and three different initial wavepacket momenta, 
r, M2 and K. Before discussing in detail, we make some 
general observations: 

• Eor the E point and M 2 point initial states, (0)pp-sc, calcu¬ 
lated using the basic ‘point-particle’ semiclassical equa¬ 
tions (11), shows similar overall qualitative features as 
the evolution of the exact (0)exact. but it generally fails 
to quantitatively reproduce the evolution. On the other 
hand, (0)wp-sc. calculated using the modified semiclassics 
of Sec. Ill (taking into account the wavepacket structure 
in momentum space) reproduces many of the prominent 
features of the (0(f))exact curve. 

Eor the K point initial state, there seems to be no notice¬ 
able agreement. 

• There is generally better agreement between the semiclas¬ 
sics and the exact evolution for the wavepacket that is 
larger in real space (cr^ = 50, top row) compared to the 
smaller wavepacket (cr^ = 0.5, bottom row). 

• Eor the larger wavepacket (cr^ = 50, top row), the E 
point initial state is almost completely in the lower band 
(W- — 0.999) and the M 2 point initial state is almost 
completely in the upper band (W+ = 0.985). Hence, us¬ 
ing only the lower band (E) or only the upper band (M 2 ) 
would give very nearly the same semiclassical curves as 
the ones shown, which are -weighted mixtures. 
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FIG. 6. (Color online) The angular velocity from exact simulations ((@)cxact), compared with point-particle semiclassics ({0)pp-sc) and with 
extended semiclassics taking momentum-space structure into account ({0)wp-sc). Initial wavepacket size: = 50 [(a)-(c)] and — 0.5 

[(d)-(f)]. Three different initial momenta are shown: ko = F [(a),(d)], ko = M2 [(b),(e)], and ko = K [(c),(f)]. For ko = F and ko = M2, 
(^)wp-sc agrees quantitatively with {0)exact, whereas (0)pp-sc shows qualitative agreement at best. For ko = K, (0)exact shows oscillations which 
are not captured by semiclassics. In panel (c), the Fourier transform of the oscillations in {0)exact is shown as inset. 


For the smaller wavepacket (cr^ = 0.5, bottom row), the 
same is true with the F point initial state {W- = 0.99), but 
the M 2 point initial state now has significant contribution 
from the lower band as well (FF+ = 0.74). This leads to 
cancellation of the weighted mixtures, so that, comparing 
panels (b) and (e), we see much smaller values of {&) for 
the smaller wavepacket. 

For the K point initial state, the contributions of the two 
bands largely cancel each other, resulting in tiny semiclas- 
sical predictions for {9{t)). 

We now discuss in more detail the larger (tr^ = 50) 
wavepacket (top row). A wavepacket with zero initial mo¬ 
mentum has support almost completely on the lower band 
(W- = 0.999), hence the relevant profile of angular velocity 
is 0-(k). At the zero momentum (F) point, both the gradient 
of the band dispersion and the Berry curvature in the lower 
band are zero, leading to a zero angular velocity. As a result, 
the wavepacket starts with zero 9. From the momentum space 
trajectory in Fig. 4(a) and the Brillouin zone profile of Fig. 5, 
one can infer that the wavepacket mostly moves through re¬ 
gions of near-zero 9. As a result, the 9 remains relatively small 
as seen in Fig. 6(a). The momentum-space shape of the packet 
plays a strong role in this case: as the trajectory lies roughly 
between positive and negative regions of 0_(k), small varia¬ 
tions of the shape can cause 9 to vary between positive and 
negative values. Accordingly, the dynamics of (0(f))exact is 
captured notably better by the extended semiclassics (0)wp-sc 
than by the point-particle approximation (0)pp.sc. 

As discussed in Sec. IIB, the wavepacket with initial mo¬ 


mentum at M 2 has support almost entirely on the upper band 
{W- = 0.015), hence the upper band characteristics are more 
relevant here. At the M 2 point, although the gradient of the 
band dispersion vanishes, the Berry curvature has a sharp peak 
[see Fig. 1(c)]. As a result, the wavepacket gains a finite angu¬ 
lar velocity almost immediately at t « 0. (Note that the semi- 
classical approximations, by construction, start with nonzero 
0 at f = 0, which is the value of 0 at the M 2 point. The 
physical or exact 0 starts at zero.) Considering the trajectory 
Fig. 4(b) and the Brillouin zone profile of Fig. 5, we note 
that the trajectory moves through regions of large 0; this is 
reflected in the larger absolute values of Fig. 6(b). The tra¬ 
jectories in momentum space intersect regions of 0+(k) < 0 
for f < 10 to explore regions 0+(k) > 0 at later times. The 
change of sign can be seen in Fig. 6(b) in all three curves. 

For the wavepacket with ko = K, the dynamics of (0) exact 
shows pronounced oscillations which preclude meaningful 
comparison with the semiclassic predictions. The oscillations 
are due to the fact that the initial state created according to 
Eq. 3 has significant weight on both lower and upper bands 
(W- = 0.457), with a relatively well-defined energy gap be¬ 
tween eigenstates occupied in the lower band and eigenstates 
occupied in the upper band. This is seen through the overlaps 
plotted in Fig. 2(a). A Fourier transform of the (0(f))exact, 
inset to Fig. 6(c), shows that the dominant frequency (peak 
around « 2.8 with width « 0.3) matches the energy differ¬ 
ence (« 2.6) between eigenstates of high overlap. Fig. 2(a) 
top right. The weighted band averages shown as semiclassi- 
cal predictions stay near zero, which could be thought of as 
the value around which the exact (0(f))exact oscillates, but it 
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FIG. 7. (Color online) Overlaps as in Fig. 2(a), for a smaller 
wavepacket, — 0.5. Comparing with the case of = 50 
[Fig. 2(a)], this smaller wavepacket has a much more spread-out dis¬ 
tribution of weights on the eigenstates of Hhu- The weight of the 
wavepackets in the lower band W- (provided with each panel), are 
all closer to 1/2 compared to the larger wavepacket of Fig. 2, where 
we had W- = 0.999, 0.015, 0.457 for these three momenta. 


is currently unclear whether this is a coincidence. It is also 
currently unclear whether a more sophisticated way of incor¬ 
porating multiple bands might allow the semiclassics to re¬ 
produce the oscillatory behavior or the average curve around 
which (0(f))exact oscillates. 

In the lower panels of Fig. 6, we have used an initial 
wavepacket with cr^ = 0.5. The Gaussian wavepacket is cen¬ 
tered at the center of a hexagon in real space, so that even with 
such a small a there are six sites equally occupied. The exact 
(0(f))exact now deviates significantly from the point-particle 
semiclassics, (0)pp.sc. The extended semiclassics, (0)wp-sc, 
continues to describe the overall behavior, for the F point and 
M 2 point initial states. The exact dynamics now shows os¬ 
cillations for all three initial momenta. This can be under¬ 
stood through the overlap distribution, shown in Fig. 7. For 
the smaller packet, the overlaps are spread out more in energy 
and also are far more biased toward more equal occupancies 
of the two bands (W_ values are closer to ^ compared to the 
corresponding values for the bigger packet). As a result, in¬ 
terference oscillations are visible also for the F point and M 2 
point initial states, panels (d) and (e). 


V. SEMICLASSICS IN A TIGHT TRAP 


In Sec. IV we showed that, as the spread of the initial 
wavepacket in real space is made bigger, the agreement be¬ 
tween the results from the semiclassical approximations and 
exact simulations improves. The idea is that, increasing size in 
real space corresponds to decreasing size in momentum space 
(as reflected also in decreasing spread in energy space); thus 
the point-particle approximation in momentum space is more 
justified. However, increasing real-space size can also lead to 
violation of semiclassics, as the semiclassic equations of mo¬ 
tion also assume sharply defined position. One effect is that 
the finite spread of the wavepacket in real space could lead 
to different potential gradients (different forces) at different 
points within the wave packet. This effect would not play a 
role for a constant gradient but can occur in a harmonic trap. 


This kind of ‘tidal’ force, makes the point particle notion less 
justified in real space. 

In order to characterize this effect, we consider the geom¬ 
etry of Fig. 8(a) with Vq = 0.02, one order of magnitude 
stronger than that in Fig. 3. The distance between the cen¬ 
ter of the trap and the center of the wavepacket is adjusted 
such that the force at the center of the wavepacket remains the 
same compared to the geometry shown in Fig. 3(a). It can be 
seen from Fig. 8(a) that, in the course of time evolution, the 
wavepacket breaks apart, spreads out and does not keep the 
notion of a well defined wavepacket as compared to the shal¬ 
lower trap [Fig. 3]. This is also reflected in the evolution of 
the wavepacket in momentum space as shown in Fig. 8(b). For 
similar time scales as those studied in previous sections, the 
wavepacket in momentum space gets distorted and diffuses 
out completely, unlike the case in Fig. 3(c) where there is still 
a notion of a well defined peak centered around some value 
of momentum. Quite surprisingly, our semiclassical calcula¬ 
tion does not seem to fail completely even in this extreme case 
since (0)wp-sc and (0)exact still agree qualitatively [see Fig. 8]. 
Even the point-particle semiclassics, (0)pp.sc, reproduces qual¬ 
itatively some of the peaks and dips of the exact curve. 

As in previous cases, for ko = K there are strong oscilla¬ 
tions due to occupancies in both bands. 


VI. SUMMARY AND DISCUSSION 

In this work, we have explored the dynamics of a gaussian 
wavepacket, with and without initial crystal momentum, on 
the Haldane honeycomb Chern-band lattice in the presence of 
external forces provided by a harmonic potential. We have fo¬ 
cused on short-time dynamics and compared to semiclassical 
treatments. Semiclassical descriptions of wavepacket dynam¬ 
ics are obviously appealing, but the range of applicability is 
poorly explored. This work can be regarded as a step toward 
obtaining detailed intuition on the regimes of validity of the 
semiclassical manner of predicting trajectories. We have for¬ 
mally treated single-particle dynamics, but our considerations 
should be applicable to non-interacting Bose condensates, and 
approximately to weakly interacting Bose condensates. 

We have found that the point-particle semiclassics re¬ 
produces many qualitative features even for wavefunctions 
that are quite small in real space, and hence extended over 
non-negligible portions of the Brillouin zone. When this 
momentum-space extent is taken into account, the agreement 
can be excellent even when the point-particle approximation 
fails. This shows that the basic idea of semiclassics, fol¬ 
lowing position and momentum simultaneously, can correctly 
embody quantum dynamics even when the point-particle ap¬ 
proach fails. Using a tight trap, we have also shown that 
this extended semiclassical approach can function reasonably 
even when the wavepacket is completely distorted or even 
torn apart (Fig. 8). While this is reassuring for the philos¬ 
ophy behind semiclassics, it does not immediately lead to a 
computationally advantageous approximation, since to obtain 
the momentum-space structure we first evolved the complete 
system in time (i.e. solved the problem computationally). 
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(a) 


(b) 






FIG. 8. (Color online) Dynamics in a tight trap, Vo = 0.02, initial position |ro — rd = 16.45, initial size: ~ 10. (a,b) Snapshots for 

k() = r. (a) Real-space occupancies (|ci(f)|^). Solid lines are equipotential contours of the trap potential, (b) Corresponding momentum 
space occupancies |ck(f)|^ in the Brillouin zone, (c) Comparison between the angular velocities calculated from full simulations ((6)exact), 
with the two types of semiclassics ((0)pp-sc and {0)wp-sc). 


However, one can envision an extended setniclassics where 
the wavepacket in momentum space is assumed to have fixed 
shape and the center moves according to the point-particle 
equations (9). 

This study is directly motivated by recent developments 
such as the experimental realization of Haldane’s honeycomb 
model in a cold-atom system [13], the experimental and the¬ 
oretical interest in the response of a localized wavepacket to 
an applied force (potential gradient) [13, 22, 27, 29, 44], and 
recent theoretical studies of dynamics in various backgrounds 
using semiclassics [29, 54, 56, 58-62]. In 2D lattices, a po¬ 
tential gradient will of course lead to Bloch oscillations in the 
gradient direction (as widely studied, e.g., in [62-68]), but 
may also induce a transverse response. This can happen even 
without Berry curvature [62], simply due to the structure of 
the energy band. In the presence of Berry curvature, the trans¬ 
verse response occurs due to a combination of the two effects. 

The present study opens up many new questions. First, we 
have focused on timescales such that the wavepacket displace¬ 
ments are of the order of the wavepacket size. It remains an 


open question to see how well semiclassics works when the 
trajectories are long compared to the wavepacket size. For ex¬ 
ample, Refs. [29, 54] have predicted a dramatic turning point 
in the trajectory of a wavepacket traveling through a honey¬ 
comb (graphene) lattice with a constant potential gradient. 
It is unclear how closely a finite-sized quantum wavepacket 
would follow such a prediction, especially when the position 
dynamics couples to internal distortion dynamics of a realis¬ 
tic wavepacket. Second, it remains an open question whether 
some version of our “extended semiclassics” can be fashioned 
into a computationally useful approximation scheme. Third, 
our treatment of multiple-band occupancies is rather prim¬ 
itive (an incoherent average), and is unable to account for 
interference oscillations. Clearly, development of multiple- 
band semiclassics for such time-evolution phenomena, per¬ 
haps along the lines of Ref. [59], is called for. Finally, 
since cold-atom experiments are more likely to track interact¬ 
ing Bose condensates or fermionic clouds rather than single¬ 
particle wavepackets, it is of interest to find out in which situ¬ 
ations the dynamics of many-fermion and many-boson clouds 
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resemble single-particle wavepacket dynamics. 
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and B turn out to be 


3 

BoCk) = 2 J 2 cos cos k • b^, 

i=l 


3 

Ba;{'k) = Ji^cosk-a^, 

i=l 


(A3a) 


(A3b) 


Appendix A: Momentum space properties of the Haldane model 


3 

f3y(k) = Ji ^ sink • a^, (A3c) 

i=l 


The semiclassical equations of motion (9) take the gradient 
of the energy dispersion and the Berry curvature in momen¬ 
tum space as input. In this Appendix, we give some details of 
these momentum space properties of the Haldane model, for 
completeness. Being realized on a lattice with a two-site ba¬ 
sis, the Haldane model Hamiltonian in momentum space can 
be written as a collection of 2 x 2 Hamiltonians /i(k): 

'Hhm = E (Al) 


with 


h(k) = a°B°{k) + a- ■ B(k) (A2) 

where 6^ j^) and Pa{b) k creates a Bloch state 

with momentum k on sublattice A{B). The ct’s are Pauli ma¬ 
trices with (T° being the identity matrix. We label the vec¬ 
tors connecting the nearest-neighbor sites as ai = (0,-1), 
^2 = ('x/S/S, 1/2), and aa = (—•\/3/2,1/2), and the vec¬ 
tors connecting the next-nearest-neighbors as bi = a 2 — aa, 
b 2 = aa — ai, and ba = ai — a 2 . With these notations, Bq 


3 

i? 2 (k) = —2J 2 sinsink • bi. (A3d) 

i=l 


The energy dispersions of the two bands are 


T;±(k)=i?o(k)±|B(k)|2, 


(A4) 


where +(—) denote the upper (lower) band. 

The eigenstates of the Hamiltonian (A2) for a given mo¬ 
mentum k can be written as 


u+,k = 


;Ck 


2 sin 


; u-,k = 




cos ^ , 


where pk and Ck are dehned via 

S.(k) 


77 k = cos 


; Ck = tan 


-1 


|B(k)|’ B^ik)- 

With these notations, the Berry curvature is given by 


= =F 


1 

— ( 

Ttt 


cos77k][5fc^Ck] 


(A5) 


(A6) 


(A7) 


where is the two-dimensional Levi-Civita and 
refers to the Berry curvature of the lower (upper) band. 
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